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We analyze simulations results of a model proposed for etching of a crystalline solid and results 
of other discrete models in the 2 + 1-dimensional Kardar-Parisi-Zhang (KPZ) class. In the steady 
states, the moments Wn of orders n = 2, 3, 4 of the heights distribution are estimated. Results 
for the etching model, the ballistic deposition (BD) model and the temperature-dependent body- 
centered restricted solid-on-solid model (BCSOS) suggest the universality of the absolute value of the 
skewness S = W3/W2'^^^ and of the value of the kurtosis Q = W/^l'W'i' — 3. The sign of the skewness 
is the same of the parameter A of the KPZ equation which represents the process in the continuum 
limit. The best numerical estimates, obtained from the etching model, are |S| = 0.26 ± 0.01 and 
Q — 0.134±0.015. For this model, the roughness exponent a = 0.383±0.008 is obtained, accounting 
for a constant correction term (intrinsic width) in the scaling of the squared interface width. This 
value is slightly below previous estimates of extensive simulations and rules out the proposal of the 
exact value a — 2/5. The conclusion is supported by results for the ballistic deposition model. 
Independent estimates of the dynamical exponent and of the growth exponent are 1.605 < z < 1.64 
and /3 = 0.229 ± 0.005, respectively, which are consistent with the relations a + z — 2 and z = a//3. 

PACS numbers: 05.40.-a, 05.50.-fq 



I. INTRODUCTION 

Surface growth processes and deposition of thin films have attracted much interest from the technological point of 
view [1-3] and motivated the proposal of continuum and discrete models for surface and interface growth, which also 
play an important role in the field of non-equilibrium Statistical Mechanics. One of the most important phcnomcno- 
logical theories is that of Kardar, Parisi and Zhang (KPZ) [4], in which the time evolution of the interface described 
by the height function h at position x and time t is given by the KPZ equation 

^ = ^V'h+^{S7hf+rj{x,t). (1) 

Here v represents a surface tension, A represents the excess velocity and ?7 is a Gaussian noise [2,4] with zero mean 
and variance (rj (x, t) r](x' , t')) = D5'^{x — x')5 {t — t'), where d is the dimension of the substrate. 

/ _2\l/2 

The interface width ^(L, t) ^ ih"^ — h ) characterizes the roughness of the interface, for growth in a substrate 

of length L (overbars denote spatial averages and angular brackets denote configurational averages). For short times, 
the interface width scales as 

e - 1^, (2) 

with P called the roughness exponent. For long times, a steady state is attained and the width saturates at 

Ut i", (3) 

with a called the roughness exponent. Eqs. (2) and (3) correspond to limits of the dynamical scaling relation of 
Family and Vicsek [5] 

e « L"/ (tL-^) , (4) 

where the dynamical exponent z = a/ (3 characterizes the crossover from the growth regime to the steady state. 
Galilean invariance gives a-\- z — 2 for KPZ in all dimensions [2,3]. The exact scaling exponents are known in d = 1, 
but no exact value was already obtained in two or more dimensions [2,3]. 

Many discrete models fall into the KPZ class, such as the restricted solid-on-solid (RSOS) model of Kim and 
Kosterlitz [6] and ballistic deposition (BD) [7]. Numerical estimates of the scaling exponents in d = 1 are consistent 
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with the exact values [6,8,9] and simulations in = 2 arc frequently used to estimate them. Most of the reported values 
of a range from a = 0.37 to a = 0.4 [6,10-13], which is confirmed by numerical solutions of the KPZ equation [14-17]. 

In 1998, assuming certain properties of the height correlation functions, Lassig obtained a quantization condition 
for the KPZ exponents which gave a = 2/5 as the only solution consistent with the range of numerical estimates [18]. 
Another consequence of his work was that the moments of the steady state height distribution. 



w^ = i^(h-hy), (5) 

would obey power-counting, i. e. they would scale as Wn ~ i"". The second moment is the squared interface width 
^ . Moreover, the validity of his assumptions requires that the steady state distribution is skewed (non-zero third 
moment), contrary to the one-dimensional case (Gaussian distribution). 

Recent numerical results of Chin and den Nijs for the RSOS and the body-centered solid-on-solid (BCSOS) models 
were consistent with a = 0.4 [19], but extensive simulations of Marinari et al for the RSOS model rule out that 
value [20]. A recent study of the KPZ equation in the mode-coupling approximation provided an estimate ^ w 1.62 
in d = 2 [21]. Despite this controversy, various works [19,20,22] confirm the power-counting property for models with 
restricted height differences (RSOS, BCSOS) and suggest that the skewness 

S = (6) 

and the kurtosis 

are universal at the steady state. According to Chin and den Nijs, these universal values do not depend on the sign 

of the coefficient A of the nonlinear term of the KPZ equation (1). 

Consequently, the value of KPZ exponents in d = 2 is still an open question, and the universality of the values 
of the skewness and the kurtosis deserves to be tested in models other than those with restricted height differences. 
The main contribution of this paper is the analysis of simulations results of a recently proposed model for etching of 
crystalline solids [23] which belongs to the KPZ class. Additional support to some of our conclusions will be provided 
by simulations results of ballistic deposition and of the temperature-dependent BCSOS model [19]. For the etching 
model, we will obtain a w 0.38 after a detailed analysis of finite-size corrections. Characteristic relaxation times 
will be calculated independently and will provide estimates of the dynamical exponent z > 1.6, while estimates of 
the growth exponent give (3 ~ 0.23. Although our error bars intercept those obtained from extensive simulations of 
the RSOS model [20], our central estimates of a are smaller and, consequently, more distant from the theoretically 
proposed value a = 0.4 [18]. On the other hand, our estimates are near those by Colaiori and Moore [21] from 
renormalization under the mode-coupling approximation. We will also show that our data for various models confirm 
the universality of the steady state skewness and kurtosis. Concerning the steady state skewness, although its absolute 
value is universal, its sign changes with the nonlinear coefficient A of the KPZ equation. 

This paper is organized as follows. In Sec. II we briefly describe the models, the simulation procedure, the methods 
to estimate Wn and the method to calculate characteristic relaxation times. In Sec. Ill, wc analyze the skewness 
and the kurtosis at the steady states. In Sec. IV, we analyze the finite-size estimates of the scaling exponents of the 
etching model, also showing some results for BD. In Sec. V we summarize our results and present our conclusions. 



II. MODELS AND SIMULATION PROCEDURE 



The model for etching of a crystalline solid of Mello et al [23] is illustrated in Fig. 1 in its growth version. The 
solids have square and simple cubic lattice structures in 1 -I- 1 and 2-1-1 dimensions, respectively. At each growth 
attempt, a column i, with current height h{i) = ho, is randomly chosen. Then its height is increased by one unit 
{h{i) ^ ho + 1) and any neighboring column whose height is smaller than ho grows until its height becomes ho- One 
time unit corresponds to growth attempts in 2 -|- 1 dimensions. In the true etching version of this model, the 
columns' heights decrease by the same quantities above. However, in this paper we will always refer to the growth 
version of Fig. 1 as "the etching model". 

In the ballistic deposition (BD) model, particles are sequentially released from randomly chosen positions above 
the substrate, follow a trajectory perpendicular to the surface and stick upon first contact with a nearest neighbor 
occupied site [2,7]. 
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In the BCSOS model defined by Chin and den Nijs [19], the substrate is a square lattice and the heights h in the 
first (second) sublattice are restricted to assume even (odd) values. Also, the nearest neighbor columns always diflFer 
in height by Ah = ±1. The energy of a given heights configuration {h} is given by E {{h}) = ■■^ jK{hi — hj)"^, 
where K is an inverse temperature parameter and the sum runs over all next nearest neighbor pairs. At each 
deposition attempt, a column c is randomly chosen and, if the constraint of the heights difference is satisfied, then 
h{c) — > h{c) + 2 with probability p = min{l,exp{—AE)), where AE is the energy change if the deposition takes 
place. The corresponding coefficient A (Eq. 1) changes sign at a critical point Kc [19,24]. 

The etching model was simulated in lattices of lengths L = 2" (n = 5 to n = 10) and lattices of lengths L = 2™ x 50 
(m = to m = 4). The maximum deposition time ranged from 10^ for the smallest lattices to 6 x 10* for the 
largest ones. During half of this time or more, the systems have undoubtedly attained their steady states. For the 
smallest lattices, 2 x 10** realizations were simulated, and nearly 300 realizations for the largest lattices {L = 800 and 
L = 1024). BD was simulated in lattices of lengths L = 2", from n = 5 (i = 32) to n = 11 (L = 2048). For the 
analysis of the skewness and of the kurtosis, it was essential to simulate this model in very large lattices, thus the 
number of realizations was relatively small: typically 10"^ realizations for L < 512, 35 for L = 1024 and 8 for L = 2048. 
The temperature-dependent BCSOS model was simulated with K = 0.25 (i = 16 to L = 128) and K = 1.0 (L = 16 
to L = 512). For L < 256, 10^ different realizations were considered, and 40 realizations for L = 512. 

The procedure to calculate average quantities, described below, followed the same lines for all models. It was 
previously used in the analysis of other growth models in 1 + 1 and 2 + 1 dimensions [9,26,27]. 

One important point is the criterion to determine the initial time tmin for estimating average quantities at the steady 
state. In this regime, the interface width ^ fluctuates around an average value instead of increasing systematically, 
which was the case in the growth and in the crossover regions. Thus, for a fixed lattice length L, the first step was 
choosing a time interval tmin < t < t„iax, with nearly constant ^, from visual inspection of the ^ x t plot (tmax was 
always the maximum simulation time). Subsequently, two tests were performed. For the first test, the interval was 
divided in 5 subintervals and the average value of ^ was calculated in each one, forming a sequence of estimates {^{i)}, 
with i = 1, . . . 5. If ^(i) < — 1) at least two times along this sequence, then the average value of ^ in the region 
tmin < i < tmax^ ^TaT^ ^ Calculated. In the second test, from the plot of log (^^™' — ^{t)) x t, for t in the crossover 
region, we obtained a rough estimate of the characteristic time r of relaxation to the steady state, as shown in Ref. 
[26]. If i„i;„ > lOr, then the interval tmin <t < t^ax was accepted as representative of the steady state. Otherwise, 
a larger value of tmin was chosen and the tests were repeated (it seldom occurred). 

In order to estimate the moments of the height distribution and their error bars, we used their average values within 
the five subintervals defined above. Final estimates of W„ are averages of these values, and error bars were obtained 
from their variances. 

The dynamical exponent was estimated using a recently proposed method [26], in which a characteristic time to, 
proportional to the relaxation time r, is calculated. For fixed L, after calculating the saturation width ^sat{L), tq is 
defined through 

^{L,To) = kUt{L), (8) 
with a constant k ^1. From relation (4), we obtain 

TO ~ L^ (9) 

Typically, the uncertainty in tq is much smaller than that of r, estimated from log {S,sat — ^ t plots - see Ref. [26]. 

We estimated To with k ranging from k = 0.5 to fc = 0.8. Although these constants are not much different, the 
values of tq for k = 0.5 and k = 0.8 typically differ by a factor 4, for fixed L. This method was already applied with 
success to calculate the dynamical exponent of other models in various universality classes [26,27]. 

III. UNIVERSALITY OF SKEWNESS AND KURTOSIS 

In Fig. 2a we show the steady state skewness S {L,t oo) of the etching model as a function of the inverse lattice 
length, which provides a good linear fit of the data, with an asymptotic estimate 5 = 0. 26 ±0.01. The absolute value 
of S agrees with results of Chin and den Nijs [19], Shim and Landau [22] and Marinari et al [20] for the RSOS model 
of Kim and Kosterlitz [6] and for BCSOS models, but those authors obtained S with negative sign. Here, the positive 
sign is related to the presence of sharp hills (see process at the left in Fig. 1) and wide valleys at the surface of the 
deposit, the opposite being observed in RSOS deposits. 



3 



Note that Eq. (1) is invariant under the transformations h ^ —h and A ^ —A, without changing the other 
parameters. This transformation changes the sign of the skewness. Consequently, we expect that the sign of S is 
related to that of A, as previously observed in the growth regimes of 1 + 1-dimensional KPZ systems [28]. In fact, in 
the true etching version of this model, with erosion leading to decreasing heights, the sign of S changes, corresponding 
to the transformation h ^ —h, X ^ —X. 

Results for the other models contribute to this discussion. In Fig. 2b we show S {L,t ^ oo) versus 1/L for the 
BD model. Note that S is negative for small lattices (typically L < 500), but for large lattices it becomes positive, 
showing that there are significant morphological differences between the steady states of small lattices and those of 
very large systems. This is the main reason to avoid extrapolating the data in Fig. 2b, even choosing extrapolation 
variables other than that 1/L (this abscissa in Fig. 2b was chosen only for illustrate the evolution of S with L). Other 
consequences of this complex finite-size behavior were previously discussed in Ref. [9]. However, it is clear from Fig. 
2b that S is asymptotically positive. 

In Figs. 2c and 2d we show the steady state skewness for the temperature-dependent BCSOS models with K — 0.25 
and K = 1.0, respectively. For K = 0.25, the extrapolation of the three last data points give S = —0.28 ± 0.015, 
reinforcing the conclusion on the universality of IS"]. Due to the small number of data points, the asymptotic correction 
term may be other than 1/L, but it would not affect that conclusion. For K = 1.0, the skewness is always positive 
and rapidly increases with L, thus no extrapolation variable of the form , with A > 0, provides a reasonable 
linear fit. However, S will certainly converge to a positive value as L — > oo, showing that the sign of the steady state 
skewness changes as K increases. 

The sign of the corresponding parameter A is obtained from the size dependence of the growth velocity. The steady 
state growth velocity, Vs (L), and the velocity in an infinitely large substrate at long times, v^a-, obey the relation [29] 

(i) = Woo - aAL-«ii , (10) 

where ay = 2 (1 — a) and a is positive. Considering a = 0.38 (see Sec. IV), we obtain ay = 1.24. In Figs. 3a-3d, 
we plotted Vg [L) versus x = l/L^-"^^ for the four models. For the etching model, the BD model and the BCSOS 
model with K = l.O (Figs. 3a,b,d), Vs decreases with x, which gives a positive A, while the opposite occurs in the 
BCSOS model with K = 0.25. It confirms that the steady state skewness of the KPZ equation has the same sign of 
the parameter A, its absolute value being universal. The same conclusion was derived from the time dependence of 
the velocity in the growth regime [29]. 

At this point, we recall that Chin and den Nijs [19] simulated the temperature-dependent BCSOS model with 
K = —0.25 and K = 0.25, obtaining S « —0.26 in both cases. They concluded that S did not depend on A, but they 
did not estimate the values of this parameter. Since A is expected to cross zero for a positive K, we conclude that 
their work failed to consider the regime of positive A, which is represented here by if = 1.0. 

Now we turn to the analysis of the kurtosis (Eq. 7). In Fig. 4a we show the steady state Q for the etching model, 
as a hmction of 1/L. The size dependence is much weaker than that of the skewness, so that extrapolation variables 
other than 1/L do not have a significant infiuence on the asymptotic estimate, Q = 0.134 ± 0.015 (the large error 
bar is mainly a consequence of the uncertainties of the finite-size data). This estimate agrees with previous ones for 
models with restricted heights differences [19,22,20], suggesting that the steady state kurtosis is also universal. In Fig. 
4b, we show the data for the BD model and for the BCSOS model with K = 0.25 and K = 1.0, which have significant 
finite-size dependence (in particular those for BD). Thus, the extrapolated values have very large error bars, but are 
still consistent with a universal value of Q. 



IV. ROUGHNESS, DYNAMICAL AND GROWTH EXPONENTS 

Our first step to estimate the roughness exponent was to calculate effective exponents defined as 

It is expected that Q!(L,i) — * oi for any value of i. 

Using different values of i in Eq. (11), we noticed that Q!(i.i) varied with L typically in the range 0.33 < Q:(L.i) < 0.38 
for 50 < L < 1024, which suggests that corrections to the scaling relation (3) are relevant. Our first proposal is to 
assume the main scaling correction as 

Ut~i"(ao + aii-^), (12) 
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where ao and ai are constants. Consequently, a{L,i) is expected to vary as 



a(i,i) «a + BL-^, (13) 

where 

^ = T^" (14) 

ln(i) ao 

Our data with 50 < L < 1024 were analyzed using four values of i in Eq. 11: i = 2, i = 2.56, i = 3.125 and i = 4 
(for noninteger i, only three or four effective exponents can be calculated). For each i, we plotted Q:(i,,i) x L^^ using 
several exponents A and did least squares fits of those plots, from which the linear correlation coefScients r(A, i) were 
obtained. Since there is no argument to predict the value of A, we adopted the condition of maximizing the coefficient 
r to extrapolate our data. In Table I, we show the exponents A which gave the largest r (best linear fits) for each i. 
The procedure is illustrated in Figs. 5a and 5b, in which we show a(L,2) x L~°-^^ and Q!(l,2) x L~^-^^, respectively, 
with the corresponding linear fits. The values of A in Figs. 5a and 5b are those which give the best fits with i = 2 
and i = 4, respectively (see Table I). 

A is expected to be independent of the particular choice of i in Eq. (11), so the differences between the estimates 
in Table I are effects of the maximization of r. Moreover, other exponents A near the values shown in Table I also 
provided reasonable linear fits of Q:(i,_j) x plots. In other words, large linear correlation coefficients were also 
obtained by considering 0.3 < A < 0.8 for different choices of i. Consequently, it is not possible to obtain a reliable 
estimate of that correction exponent. On the other hand, the asymptotic a obtained from the same fits fluctuate 
within a narrow range. Accounting for the error bars of the data, we obtained a = 0.385 ± 0.01 for « = 2 and 
a = 0.382 ± 0.01 for i = 4. 

We also checked the effect of considering a fixed correction exponent A = 0.55 (which is near the values in Table I) 
for all values of i. This procedure will be called fixed A method. Least squares fits of the a(L.i) x i^O-S'' plots were 
performed, providing the slopes B which are shown in Table I. From Eq. (14), it is expected that ^^_'"a^ ~ ^ ~ const 
(see also Eqs. 12 and 13), thus we also showed in Table I the corresponding estimates of ^^_'"a) • This quantity 
fluctuates with i, indicating that there is no systematic trend in our results due to choice of different values of i for 
calculating effective exponents. Also notice that the estimates of a from the fixed A method are in same range of the 
those obtained from the maximization of correlation coefficients. 

Since the range of lattice lengths considered here is not very large and the correction exponents A estimated above 
are relatively small, we tried to improve our analysis with a different assumption for the scaling corrections. Contrary 
to the previous procedure, now we will consider a well defined form for the main scaling correction, which is an 
additional constant term ^| in the dynamic scaling relation (4) for the squared interface width: 

e{L,t)=ei+L'''9{tL-^), (15) 

where 5 is a scaling function. ^7 is called intrinsic width and represents contributions of small length scales fluctuations, 

typical of models with large local heights differences [30-32] such as the etching model and the BD model. From 
Eqs. (15) and (12), the assumption of the intrinsic width as the most relevant subleading correction corresponds to 
a (fixed) correction exponent A = 2a w 0.8. It is slightly larger than the typical values obtained in the previous 
analysis. 

Effective exponents a^j^^ which cancel the contribution of ^| are defined as 



^(/) ^ 1 In [e^at (2L) - esat (L)] I {Hat (L) ' iL/2)] .^g. 
"^2 In 2 ■ ^ ' 



In Fig. 6a we show a'^j^' versus 1/L for the etching model and a least squares fit of these data, which provides 

a = 0.383 asymptotically. Here, the variable 1/L in the abscissa was chosen only to illustrate the behavior of the 
data. It represents a second correction term for the ^ scaling, which is still more difficult to measure than the first 
correction term. Thus, we also tested other variables in the form to extrapolate a^P, with 0.5 < Ai < 2 (Ai = 1 

was used in Fig. 6a). The corresponding linear fits give 0.380 < a < 0.387. The small range of the asymptotic a is 
a consequence of the slow variation of a^P with L (2% from L = 100 to L = 500), as shown in Figs. 6a. Accounting 
for the error bars of the data, our final estimate is a = 0.383 ± 0.008. 

The forms of finite-size corrections analyzed above, Eqs. (13) and (15), cannot be rigorously justified, but are based 
on heuristic arguments. See e. g. Ref. [2] and references therein. Certainly, the fact that increases slowly with 
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L is a support to the assumption that the intrinsic width is the most relevant correction term to Anyway, using 
different assumptions on the scahng corrections was essential to confirm the reliability of the above estimate of the 
roughness exponent. 

The universality of the values of skcwness and kurtosis imply that W3 and W4 may also be used to estimate a. The 
effective exponents obtained from the third moment have very large error bars, but those obtained from the fourth 
moment behave similarly to the ones obtained from the interface width. The analogous of exponents ct(L,'2) and q;(l,4) 
(Eq. 11) calculated with W4 also converge to the range 0.38 < a < 0.385 with strong corrections to scaling. Effective 
exponents which cancel the contribution of a constant additive term in the dynamic scaling relation for W4 (analogous 
to a^j^^) are defined as 



(/,4) _ 1 In [W4 (2L, t ^ 00) -W4{L,t^ 00)] / [Wj {L, t ^ 00) - Wj (L/2, t ^ (X))j 
-4 ln2 • 



They are plotted in Fig. 6b as a function of 1/L. The asymptotic estimate, obtained with the procedure described 
above, is a = 0.379 ± 0.012. It is in good agreement with the estimate from the interface width and also excludes 
a = 0.4. 

The same analysis was also performed with the BD model, as shown in Fig. 7, with data for L < 512. Although 

the results are less accurate than those for the etching model, they also suggest that a < 0.4 asymptotically. 
In order to estimate the dynamical exponent z, we calculated effective exponents zl defined as 

_ ln[ro(2L)/ro(L/2)] 
= h^4 ' 

using the characteristic times tq defined in Sec. II (Eqs. 8 and 9). In Figs. 8a and 8b we show zl versus 1/L for the 
etching model, obtained using k = 0.6 and k = 0.8 to calculate tq, respectively. Here, the abscissa 1/L is also chosen 
to illustrate the evolution of z^, but not to perform extrapolations of the data. Although the error bar of the data for 
L = 800 is relatively large, those plots indicate that z > 1.6. Considering the trend for large L and different values 
of k, we estimate 1.605 < z < 1.64. These values are consistent with the above estimates of a and the exact relation 
a + z = 2. 

Finally, we estimated the growth exponent (Eq. 2) of the etching model with the same procedure previously applied 
with success to BD and to the Das Sarma and Tamborenea model in 1 + 1 dimensions [9,27]. The growth region 
for each L begins at to = 50 and ends at the maximum time Tmax such that the linear correlation coefficient of 
the data in the range to < t < Tmax exceeds a fixed value r-min [9]. Here, rmin = 0.99995 and rmin = 0.9999 are 
considered. Effective exponents /3l are defined as the slopes of the linear fits of logPT x logt plots using all data 
in the above-defined growth regions. In Fig. 9 we plot fS^ versus 1/L for the etching model using rmin = 0.99995 
and rmin = 0.9999. The error bars of those effective exponents is very small. Wc also show in Fig. 9 the linear fits 
of the /3l X 1/L data for each rmin- Other variables in the form L~^ were used to extrapolate the Pl data, giving 
asymptotic estimates /3 0.229 ± 0.005. Within error bars, it agrees with the value a/z = 0.234 ± 0.009 obtained 
from the above estimates of a and z. 



V. CONCLUSION 



We studied four 2 + 1 dimensional discrete growth models in the KPZ class, determining critical exponents and 
steady state values of the skewness S and the kurtosis Q. Accurate estimates of the scaling exponents were obtained 
for the etching model proposed by Mello et al [23]: a = 0.383 ±0.008, 1.605 < z< 1.64, P = 0.229 ±0.005. Results for 
the ballistic deposition model also indicate that a < 0.4. The presence of the intrinsic width as the main correction 
to the interface width scaling was considered to extrapolate the simulations data. We also obtain the estimates 
S = 0.26 ± 0.01 and Q — 0.134 ± 0.015 in the steady state regime of the etching model. Results for the BD model and 
of the BCSOS model, together with previous results for the RSOS model, suggest that the absolute value of S and 
the value of Q are universal, the sign of the skewness being the same of the parameter A of the corresponding KPZ 
equation. 

The above estimate intercepts the error bar of the roughness exponent of the RSOS model given by Marinari and 
co-workers [20], a = 0.393 ± 0.003. However, our central estimate is significantly lower than theirs, and the result for 
BD confirm this trend. On the other hand, our estimate is very near that by Colaiori and Moore, a « 0.38, from 
renormalization methods [21]. Additional support to our conclusions was given by the independent calculation of 
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exponents z and contrary to recent simulation works on these lines [19,20], which were limited to the calculation 
of exponent a. 

We believe that much more accurate estimates of a are difficult to be achieved with numerical simulations of this 
type of lattice model. However, wc consider that this work provides a significant amount of numerical results indicating 
that the theoretically proposed value a = 0.4 [18] is not valid, within the limits of the assumptions made about the 
form of finite-size scaling corrections. This result and the evidence of universality of the values of the skewness and 
the kurtosis may motivate further analytical (maybe also numerical) studies of the KPZ theory in 2 + 1 dimensions. 
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TABLE I. For each i in Eq. (11), are given: the values of the correction exponent A which provided the largest correlation 
coefficients r of fits of Q!(L,i) x and the corresponding asymptotic a, the values of the slope B of the fits with fixed A = 0.55 
and the corresponding values of asymptotic a and of Blni/{1 — i^). 



(Eq. 11) 
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a 


B\iii/{l-i'^) 




(largest r) 
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FIG. 1. The growth rules of the model of Mello et al (original growth version). The arrows indicate the randomly chosen 
columns and the new particles are shown in gray. 

FIG. 2. Steady state skewness of various models versus inverse lattice length: (a) etching model (in the growth version); (b) 
ballistic deposition; (c) BCSOS model with K = 0.25; (d) BCSOS model with K = 1.0. The solid lines in (a) and (c) are least 
squares fits of the data. 

FIG. 3. Steady state velocities versus 1/L"ll, with a\\ = 1.24, for: (a) etching model; (b) ballistic deposition model; (c) 
BCSOS model with K = 0.25; (d) BCSOS model with K = 1.0. Error bars are smaller than the size of the data points. 



FIG. 4. Steady state kurtosis of various models versus inverse lattice length: (a) etching model; (b) ballistic deposition model 
(squares), BCSOS model with K = 0.25 (triangles) and BCSOS model with K = l.Q (crosses). Error bars in (b) (not shown) 
are smaller than the size of the data points, except for BD in the largest lattice (L = 1024). 



FIG. 5. Effective exponents a(L,2) and a(L,4) for the etching model versus L ^, with the exponents A that give the best 
linear fits for i = 2 and i = 4, respectively. 

FIG. 6. Effective exponents a^' (from interface width) and a^'*' (from the fourth moment of the height distribution) for 
the etching model versus inverse lattice length. 



FIG. 7. Effective exponents a}^ ' for the ballistic deposition model versus inverse lattice length. 



FIG. 8. Effective exponents zl for the etching model versus inverse lattice length. Relaxation times To were obtained with: 
(a) k = 0.6 and (b) k = 0.8 in Eq. (8). 

FIG. 9. Effective exponents Pl for the etching model versus inverse lattice length, obtained with minimum correlation 
coefficients rmin = 0.9999 (squares) and rmin = 0.99995 (triangles). Solid lines are least squares fits of each set of data 
(200 <L< 1024). 
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